Origin of rebounds with a restitution coefficient larger than unity in nanocluster 

collisions 
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We numerically investigate the mechanism of super rebounds for head-on collisions between nan- 
oclusters in which the restitution coefficient is larger than unity, ft is confirmed that the temperature 
and the entropy of the nanocluters decrease after the super rebounds by our molecular dynamics 
simulations, ft is also found that the initial metastable structure plays a key role for the emergence 
of the super rebounds. 
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I. INTRODUCTION 

Nanoclusters are technologically important for the con- 
struction of nanodevices. Because the size of nanoclus- 
ters is mesoscopic, thermodynamic properties of such ma- 
terials are still not well understood[l|, though the meth- 
ods to make nanoclusters such as adiabatic expansion 
through a nozzle and a laser ablation technique are well 
established. Q 

Dynamics of nanoclusters are extensively investigated 
from both scientific and technological interest. There are 
many numerical studies on cluster-cluster and cluster- 
surface collisions based on the molecular dynamics simu- 
lation. We observe variety of rebound processes for 
such systems caused by the competition between the at- 
tractive interaction and the repulsive interaction of two 
colliding bodies. Binary collisions of identical clusters 
cause coalescence, scattering, and fragmentation depend- 
ing on the cluster size and the impact energy. [TJ [l2| 
On the other hand, cluster-surface collisions induce soft 
landing, embedding, and fragmentation. 

0, EE 13 The 

attractive interaction plays crucially important roles in 
collisions of nanoclusters, so that the modeling of the co- 
hesive collisions in various scale are actively discussed in 
these days, [ill-ill . 

However, the attractive interaction can be reduced by 
some combinations of the two interacting objects and the 
relative configuration of colliding molecules. [22| Awasthi 
et at introduced a modified Lennard- Jones model to 
simulate the rebound process of a Bi cluster onto a SiC>2 
surface [10], in which they introduced a cohesive param- 
eter to reduce the attractive interaction between differ- 
ent atoms on the surface. There exists a correspond- 
ing experiment on the rebound process of Bi nanoclus- 
ters on a Si surface with the aid of a V-shaped template 
etched on a silicon wafer. [3] This suggests that the mod- 
ified Lennard- Jones clusters can approximately describe 
collisions of real clusters. Similarly, recent papers have 
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reported that surface-passivated Si nanoclusters exhibit 
elastic rebounds on Si surface due to the reduction of the 
attractive interaction between the surfaces. p| [23|, [24| In 
particular, Saitoh et al. confirmed that the behavior in 
collisions of modified Lennard- Jones clusters is similar to 
that of H-passivated Si clusters from their simulation. [3] 
These results also support that the modified Lennard- 
Jones model can be regarded as a simplified model of 
nanoclusters. 

In general, thermal fluctuation also plays an important 
role for small systems such as nanoclusters. Indeed, the 
present authors carried out the molecular dynamics simu- 
lation of colliding thermally activated modified Lennard- 
Jones clusters to investigate the impact phenomena and 
found the existence of the "super-rebound" in which the 
restitution coefficient is larger than unity. @ Recently, it 
has been reported that such a rebound can be observed 
in a molecular dynamics simulation of a collision of Cu 
nanoparticles on a rigid wall.[5| In addition, another re- 
search group has reported that large recovery strain after 
loading (called super-elasticity) can be found in a small 
system. [25} Although there is a possibility that such an 
anomalous mechanical property of small systems is also 
concerned with the emergence of "super rebounds" , the 
underlying mechanism for the "super rebounds" is still 
unclear. 

In macroscopic systems, the restitution coefficient 
larger than unity can be observed in oblique collisions 
between a hard sphere and a soft elastic plate [2Q, [27| , 
which differs from the super rebound. Indeed, the resti- 
tution coefficient can easily become larger than unity if 
the rebound direction is changed in oblique collisions. As 
another example, a recent study reported that the large- 
scale magnetized plasmoids can gain the increased kinetic 
energy after collision in the heliosphere.[28j On the con- 
trary, in microscopic systems such as nanoclusters, the 
true "super rebounds" can be observed in normal col- 
lisions, which implies that a part of thermal energy is 
converted to macroscopic degrees of freedom. This may 
imply a possibility to make a nanoscale object which ex- 
tracts work from its internal energy. Although the av- 
eraged behavior of the restitution coefficient against the 
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relative impact speed can be approximately understood 
by the macroscopic theory of cohesive collisions, [f|[l9| the 
mechanism of curious energy transfer against intuition of 
the thermodynamics is important to both fundamental 
physics and energy technology. 

The aim of this paper is to study the mechanism of the 
super rebounds from the viewpoint of the energy transfer 
between the microscopic and the macroscopic degrees of 
freedom. For this purpose, we investigate the characteris- 
tics of thermodynamic functions such as temperature and 
entropy in the super rebounds based on the molecular dy- 
namics simulation. In addition, we also investigate the 
change in structures of the colliding nanoclusters charac- 
terizing super rebounds by the introduction of Steinhart's 
order parameter 29] and some related geometrical order 
parameters. 

The organization of this paper is as follows. In the 
next section, we introduce our nanocluster model. In 
Sec. Ill we show the results of our molecular dynamics 
simulations. In Sec. IV we discuss the reason why super 
rebounds can be observed in our simulation. Section V is 
devoted to the summary of this work. In Appendices A 
and B, we explain the method of calculating entropy in 
our simulation and how to calculate Steinhardt's order 
parameters, respectively. 
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where is the distance between two arbitrary atoms 
i and j in this system, a and e are the core diameter 
and the scale of energy, respectively. Here, c is the cohe- 
sive parameter which changes the magnitude of cohesion 
between atoms. In our simulation, we adopt c = 0.2 
between the atom i on the lower surface of C p and the 
atom j on the upper surface of C while c = 1 between 
all the atoms in each cluster. Here the definition of the 
surface and the bulk atoms is as follows. Each cluster 
consists of 9 layers of atoms. We define the atoms out- 
side the concentric inscribed sphere as the surface atoms 
while the other atoms as the bulk atoms. We assume 
that the surface atoms are different from the bulk atoms, 
which may cause the weak cohesive interaction between 
the clusters. [13 HI] 

Let us introduce the volume fraction ip which is defined 
by 
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where p is the number density. Figure [2] is the relation- 



II. MODEL 

Let us introduce our model. Our model consists of 
two identical nanoclusters, each of which contains 236 
"atoms" (Fig. [I}. The clusters have facets due to the 
small number of "atoms" . To construct one cluster, we 
first make the face-centered cubic crystal of 9 x 9 x 9 
layers of atoms. Next, we cut out a spherical cluster 
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FIG. 1: (Color online) Snapshot of our model after equilibra- 
tion to T — 0.02e. Each of them contains 236 "atoms" which 
are bound together by the modified Lennard- Jones potential. 

from the cube. In FiglTJ we call the upper projectile and 
the lower target clusters C p and C', respectively. All the 
atoms in one cluster are bound together by the modified 
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FIG. 2: Relationship between volume fraction of atoms and 
potential energy of one cluster. 

ship between <p and the initial potential energy of one 
cluster. To construct a cluster, we adopt the value for ip 
as ip = 0.56 in which the potential energy is minimum, 
while we adopted <p — 0.21 in our previous work. Q 

The procedure of our simulation is as follows. As the 
initial condition of the simulation, the centers of mass of 
C p and C* are placed along the z— axis with the separa- 
tion 4i? between the centers of mass of C p and C* , where 
R is the radius of the cluster (Fig. [IJ. The two clusters 
are placed in mirror symmetric positions with respect to 
z = so that a facet of one cluster is placed face-to- 
face with that of another cluster. We have checked that 
the relative rotational orientation of the clusters around 
z-axis little affect the relationship between the restitu- 
tion coefficient and the colliding speed in our previous 
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work. Q The initial velocities of the atoms in both C p 
and C* satisfies the Maxwcll-Boltzmann distribution with 
the initial temperature T. The sample average is taken 
over different sets of initial velocities. 

Wc numerically solve the equation of motion of each 
atom i described by 
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dt 2 



(3) 



where M and x$ are the mass of an atom and the posi- 
tion of the j-th atom, respectively. is the modified 
Lennard-Jones force which is calculated from Eq. (flj as 
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Numerical integration of the equation of motion for each 
atom is carried out by the second order symplectic in- 
tegrator with the time step dt — 1.0 x 10~ 2 a j \j e j 'M . 
To reduce computational costs, we introduce the cut-off 
length r c of the Lennard-Jones interaction as r c = 2.5a. 
The rate of energy conservation, \E(t) — E 1 / \E | , is kept 
within 10~ 5 , where E is the initial energy of the system 
and E(t) is the energy at time t. 

We equilibrate the clusters by using the velocity scal- 
ing method [3(J Hl| as the thermostat in the initial 2000 
simulation steps. Here, we introduce the kinetic temper- 
ature, 



t = — y -M( Vj 



v c ) 2 , 



(5) 



where N and v c are the number of atom and the velocity 
of the center of mass of one cluster, respectively. It has 
been confirmed that the temperature of the system con- 
verges on a desired temperature during the equilibration. 

After the equilibration we give the translational ve- 
locity to C p and C* to make them collide against each 
other without the thermostat. We give the translational 
speed with the acceleration g — 0.02e/(aM). The typi- 
cal value of the relative impact speed in our simulation 
is V = 0.1y/e/M, which is slightly less than the thermal 
velocity of the system defined by Vth = y/T/M when 
T = 0.02e. 



III. SIMULATION RESULTS 
A. Macroscopic Properties 

Here we show the macroscopic properties of the collid- 
ing nanoclusters in our simulation. To characterize the 
rebound behavior of macroscopic inelastic collisions, we 
calculate the restitution coefficient e defined by 



V z (t*) 
V z (0) ' 



where V z (i*) is the z component of the relative trans- 
lational speed of C p to C*, and t* is t he sc aled time of 
separation of the clusters, i* =t*/(a/ y/e/M). We define 
the time of separation t* by the time when the relative 
rebound speed becomes stable after the collision. 
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FIG. 3: (Color online) Histogram of e for several tempera- 
tures of clusters with V z = 0.2y/e/M. 

Figure [3] shows the histograms of e for initial temper- 
atures T — 0.01e,0.02e, and 0.04e, respectively. Each 
histogram is constructed from 1000 samples with the ini- 
tial speed of V z = 0.2-^/e/M. We find 12 samples of the 
super rebound at T — 0.04e while all samples are the 
ordinary rebounds at T = O.Ole. At T = 0.02e we found 
that only 1 sample is the super rebound. On the other 
hand, if we adopt V z = 0.3-^/e/M, most of the rebounds 
are ordinary. We can find only two samples of super re- 
bounds even at T = 0.04e. Thus, the super rebounds can 
be observed only if the colliding speed is lower or almost 
equal to the thermal speed. 



B. Thermodynamic Properties 

Here we show the time evolution of kinetic temperature 
defined by Eq. ([5]) . Figure |4] shows typical examples of 
the changes in temperature in an ordinary rebound and 
a super rebound, t is the scaled time, t = tj (a/y/e/M). 
A slight discrepancy in temperature between C p and C* 
is observed after the equilibration in Fig. @] (a). The 
temperature increases after the collision around t = 9 for 
an ordinary rebound (see Fig. 0Ja)), where a part of the 
translational energy is converted to the internal energy. 
On the other hand, the temperature decreases in C* while 
the temperature of C p increases after the collision for a 
super rebound (Fig. Sfb)). 

Next let us investigate the entropy change AS defined 
by Eq. (|A9[) in Appendix [Al Figure [SJa) plots the time 
evolution of AS* for C p and C* in a super rebound, which 
shows that the entropy increases in C* while it decreases 
in C p during collision for t < 15.8. On the other hand, 
Figurejn^b) shows that the discrepancy of entropy change 
remains finite after the collision for i > 15.8. In addition, 
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FIG. 4: (Color online) Time evolution of temperature for (a) 
an ordinary rebound and (b) a super rebound. 

it is remarkable that the entropy of C* decreases corre- 
sponding to the decrease of temperature. 

So far, we have reported the results of our simulation 
on thermodynamic quantities. What we confirm is that 
the behavior of such quantities in super rebounds is seem- 
ingly in contrast to what is expected from the second law 
of thermodynamics. However, we should note that the 
second law is only strictly valid after ensemble average is 
taken over. In this sense, our result does not violate the 
second law of thermodynamics. 



C. Structural Change during Collision 

Let us investigate the structural change of the clus- 
ters induced by collisions. To characterize the structural 
change, we first introduce a local bond order parameter 
known as time-averaged Steinhart's bond order parame- 
ter Qi (see Eq. (|B3p for its explicit definition). (29l.l32l.l33j 
We note that Steinhardt's bond order parameters are 
used to characterize the structural change in nanoclusters 
associated with collisions [13] and melting [35|, [36[ as well 
as equilibrium structures of crystalline solids. According 
to the definition of the bond order parameter summa- 
rized in Appendix [B] we calculate Qi{i) and Qe(i) of 
each atom i before and after collisions, respectively. For 
the time average in calculation of Qi(i) and Qe(i), we 
use % = 5.5a/iJe/M and t a — 20a/^/e/M for before 



FIG. 5: Time evolution of AS of a super rebound in (a) 
i < 15.8 and (b) i > 15.8. Chain and solid lines show AS 
in C p and C*, respectively. At t = 30, AS of C* becomes 
-0.243. 

TABLE I: Numbers of FCC "atoms" for ordinary rebound. 

numbers of FCC "atoms" Q4, Qa 
before collision 649 0.19105 0.56897 

after collision 574 0.1908 0.5689 



and after collisions, respectively, with the time interval 
r a = l.bcr/y/ e/M, where t a and tj respectively corre- 
spond to t in Eq. (|B3|) . 

On Q^-Qe plane, the highest peak can be found around 
(Q4, Q&) = (0.190, 0.574) which is characteristic for FCC 
crystal structure. We cannot find visible shift of the peak 
position on Q^-Qe plane during the collision in both su- 
per and ordinary rebounds. However, we find that the 
number of "atoms" at the highest peak decreases after 
the collision for ordinary rebounds, which means that 
the number of "atoms" in FCC bond order decreases af- 
ter the collision (see Table QJ. On the contrary, for super 
rebounds, the number of "atoms" at the highest peak 
increases after the collision, which means that the FCC 
bond order becomes intensive after the collision. Thus, 
the super rebounds are characterized by the increase of 
the number of "atoms" with FCC bond order during col- 
lisions. 

To investigate the details of the structural change, we 
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TABLE II: Numbers of FCC "atoms" for super rebound. 

numbers of FCC "atoms" Q4 Qe 



before collision 
after collision 



690 
718 



0.1890 0.5735 
0.1911 0.5692 



also calculate the distribution of | (fori W 1 2 defined by Eq. 
(Bl). To construct these distributions, we use the parti- 
cle data at the final instant of the initial thermalization 
by the velocity scaling method. Because qi m is equal to 
(— l) m g ; * m , we only investigate the case of m > 0. 

We present the histogram of |f?64(*) | 2 in Fig. [5] The 
solid curve represents the frequency distribution for or- 
dinary rebounds while the broken curve is that for super 
rebounds. One can find the slight discrepancy between 
those two distributions. We should mention that the dis- 
crepancy is invisible for the other m, though we do not 
present the corresponding figures. 




0.12 



FIG. 6: Histogram of \q 64 \ 2 . 

To quantify the discrepancy between those two distri- 
butions, we introduce Xm defined by 



V (OrJJ) ~ S m (j)) 2 

*f o m (j) + s m {j) 



E^iW. (7) 



where O m (j) and S m (j) are the frequencies of IgwC*)! 2 
in the j-th bin for ordinary and super rebounds, respec- 
tively. We calculate Xm f° r each value of m to investi- 
gate the relationship between m and Xm- Figure [7] shows 
the relationship between m and Xmi which characterizes 
the difference between super rebounds and ordinary re- 
bounds. This figure shows that the discrepancy between 
the distributions of |(?64(*)| 2 is remarkably large compared 
to other values of m. In addition, we show the relation- 
ship between 74 (j) introduced in Eq.© and | Qe4 (*) 1 2 in 
Fig. [SJ where a remarkable difference can be found in the 
range of \q 6 4(i)\ 2 < 0.02. 
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FIG. 7: Relationship between m and Xm which is defined by 
Eq. ©. 
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FIG. 8: 7I as function of |g64(i) 



IV. DISCUSSION 

Let us discuss our results. In the previous section, 
we have demonstrated that the structural difference of 
clusters can be found between the super and the ordinary 
rebounds. Here we investigate the potential energy of 
the characteristic local structures included in the clusters 
which induce the super rebounds. 

Focusing on the noticeable peaks in Fig. |8j we select 
the "atoms" with the order in the range, 1.1643 x 1CP 2 < 
1964(0 1 2 < 1.7425 x 10~ 2 , from 20 clusters which induce 
the super rebounds. The reason why we choose the range 
of I (764 (*) 1 2 is as follows. The characteristic peaks in Fig. 
[8jare characterized by \q 6 4(i)\ 2 = (1.4534±0.2891) x 1CT 2 , 
so that we pick up "atoms" in the above range. For 
each of the selected "atoms" , we define the neighboring 
"atoms" which are located within the distance of 1.6a. 
We define a local structure by the collection of the se- 
lected atom and its neighboring "atoms" . Figure [H] is a 
snapshot of a cluster including two local structures, each 
of which is centered by the atom which has the order (red 
or medium grey). This figure shows that the local struc- 
tures are located on the surface of the cluster. Moreover, 
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the local structures seem to be unstable due to the less- 
ordered orientation of "atoms", which may lead to the 
high potential energy. 




FIG. 9: (Color online) Cluster which induces super rebounds. 
Red (medium grey) particles show the "atoms" with the order 
1.1643 x 10" 2 < |q 6 4(i)| 2 < 1.7425 x 1(T 2 while blue (dark 
grey) particles are neighboring "atoms". 

Next, we calculate the averaged potential energy of 
those local structures. For the purpose, we calculate the 
potential energy of each atom in the cluster with the 
cutoff 3.5cx and sum up the potential energy over all the 
"atoms" in a local structure to obtain its potential en- 
ergy. Among 20 clusters which induce super rebounds, 
we find 372 "atoms" to satisfy 1.1643 x 1CT 2 < \q ei \ 2 < 
1.7425 x 10~ 2 , by which we calculate the averaged poten- 
tial energy of one local structure. From our calculation, 
the averaged potential energy per one local structure be- 
comes -97.827e. 

For comparison, we focus on the potential energy of 
the local structure which belongs to the highest peak 
in Fig. |6] because we find same number of the local 
structures both in super rebounds and in ordinary re- 
bounds. The local structure in this case is defined by 
the collection of the "atoms" with the order around the 
peak, 2.8 x 10~ 2 < |g 64 | 2 < 3.0 x 10~ 2 , and its neighbor- 
ing "atoms" in the clusters which induce the super and 
the ordinary rebounds. From the calculation with 482 
local structures in the clusters which induces super re- 
bounds, we obtain the averaged potential energy per one 
local structure as — 104.814e. In the case of the clusters 
which induce ordinary rebounds, the averaged potential 
energy becomes — 104.202e from the calculation of 512 
local structures. 

From these results, the clusters for the super rebound 
has more metastable local structures with high potential 
energy rather than the clusters for the ordinary rebound. 
Thus, we conjecture that a part of the high potential en- 
ergy is transferred to the macroscopic degrees of freedom 
during the super-rebound collision. Moreover, we can 
also conjecture that the decreases of temperature and 
entropy after collision are connected with the decrease 
of the potential energy. In this paper, we focus on the 
relationship between the local structure and the super 
rebound. It is also important to understand the role of 
macroscopic structural change. The macroscopic defor- 
mation of nanoclusters after the collision will be reported 
elsewhere. [3 7j 



V. CONCLUSION 

In this paper, we have performed the molecular dynam- 
ics simulation of colliding nanoclusters to study the ori- 
gin of super rebounds. Through the investigation of the 
thermodynamic properties, we have found that the de- 
creases of the temperature and the entropy are observed 
in one of the binary clusters when the super rebound is 
observed. This may be attributed to the biased distribu- 
tion of the local bond order parameters after the collision 
which is caused by the initial metastable configuration of 
"atoms" . 

Through the investigation of the local bond order pa- 
rameters, we have found the discrepancy of the distri- 
butions of the order | (?e4 (*) 1 2 between the super and the 
ordinary rebounds. The averaged potential energy of the 
local structure which belongs to the characteristic peak 
in the |<7e4 (*) 1 2 distribution is larger than the typical lo- 
cal structures which are abundant in our cluster. Thus, 
the energy transfer from the local structure to the macro- 
scopic degree of freedom accompanied by the structural 
change during collision play an important role in super 
rebounds. Moreover, the change of entropy associated 
with the change in the internal state of order after the 
collision is the key point for understanding the emergence 
of the super rebounds. Further investigation of the re- 
lationship between the local structural change and the 
entropy decrease will be one of our future tasks. 
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Appendix A: Calculation of entropy 

In this appendix, we explain the method to calculate 
the entropy in our system. Let us assume that the col- 
lision of C p onto C* does the work 5W and changes the 
internal energy of C*. Because each cluster is an isolated 
system, we assume that the change in heat 5Q of C* is 
determined by the first law of thermodynamics, 

5Q = dE + 5W, (Al) 

where dE is the change in the internal energy. 
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The internal energy E is calculated as 

i 

2m 



E 



(A2) 



where Pi, m are the relative momentum of an atom in 
C* against the center of mass and the mass of an atom, 
respectively. U (ry ) is the potential energy between the 
"atoms" i and j in C*. 




FIG. 10: Interaction between an atom in C and an atom in 
C p . 

The amount of work 5W(t) acting on C* is defined as 
follows. Let us define the work done by the atom j of C p 
to the atom i of C* (Fig. ITUl) . Assuming that the dis- 
placements of the "atoms" are caused by the interactive 
force Fij , the amount of work 6 Wi j (t) done by the atom 
j to the atom i during an interval dt may be written as 



[r tf (t + dt) -rij-(t)], 



(A3) 



where Tij(t) = r»(t) — Tj(t). Thus, the work SW(t) done 
by C p to C* may be expressed as 



(A4) 



where we adopt the simple trapezoidal rule for the eval- 
uation of the integral. 

TABLE III: Bond orientational order parameters for fee, bec, 
and hep crystals[33[. 



<1(> 



fee 0.191 0.575 
bec 0.036 0.511 
hep 0.097 0.485 



Appendix B: Calculation of order parameter 



In this appendix, we introduce Steinhardt's order 
parameter to characterize the structure of our model 
[29l I32I l33j |. Let us assume that an atom i in a crys- 
talline structure is surrounded by Nb(i) "atoms" within 
the cutoff length r c . That is, the distance ry between 
the i-th atom and one of the neighboring "atoms" j is 
less than the cutoff length r c . In our analysis, we adopt 
r c = l.6<7. 

First, for each atom, we calculate the average of spher- 
ical harmonics which depends on ry as follows: 



N„(i) 



qim(i) 



N b (i) 



(Bl) 



Next, we introduce qi(i) according to the following defi- 
nition: 



Using these quantities, we define the entropy difference 
AS(t) = S(t) - 5(0) as 



AS(t) 



SQ(t ) 
T(t') 



(A5) 



Here we use the kinetic temperature Eq.([5} as T. To 
calculate AS(t) by Eci. ([A5[) . we first smooth the temper- 
ature and the heat as 



f(t) 
6Q(t) 



t+At 

t-At 
t+At 

t-At 



T(t )dt 



5Q(t )dt 



(A6) 

(A7) 
(A8) 



with At = .565cr / \j e / 'm because those data change 
rapidly against time. By using those smoothed parame- 
ters, the entropy change of C* may be rewritten as 



AS(t) 



6Q(t ) 

f(0 



(A9) 



\ 



An 
21 + 1 



\qim(i)\ 



(B2) 



m——l 



Instead of qi(i) 7 we use the time-averaged value defined 
by 



Qi(i) = - 



'a Jt 



ki+r a 



qi(i)dt 



(B3) 



with the time interval r a . 

Each crystalline structure, such as a body-centered cu- 
bic (bec) and a face-centered cubic (fee), has a charac- 
teristic distributions of those parameters. For example, 
the distributions of (74 and in a fee structure shows two 
peaks while the distributions in other structure has a sin- 
gle peak in a Lennard- Jones system. [33j The mean values 
of those distributions are summarized in Table Mil 32] 
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